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1.  Introduction.  We  consider  here  some  results  on  the  existence  and  the 
asymptotic  behavior  of  solutions  of  the  singularly  perturbed  nonlinear 
problems 

(1.1)  eAu  = p,A(x,u)*  vu  + h(x,u) 
and 

N 2 

(1.2)  eAu  = \j,  Z b (x,u)u  + h(x,u) 

j=l  J J 

N 

for  x in  a bounded  region  in  B (N  > 2)  which  satisfy  along  the  boundary 
of  n either  Dirichlet  - or  Robin  - type  conditions.  The  parameters  e and  y, 
are  assumed  throughout  to  be  small  and  positive  and  to  tend  to  zero  in  an 
interrelated  fashion.  We  will  determine  precisely  how  the  relative  sizes  of 
e and  |jb  influence  the  behavior  of  solutions  of  (1.1)  and  (1.2),  especially 
near  the  boundary  of  Q . At  the  same  time  we  will  give  estimates  on  these 
solutions  in  regions  where  up  to  now  complete  information  has  been  lacking. 

Such  regions  often  have  corners  and/or  boundaries  which  coincide  with  solutions 
of  certain  lower-order  equations  (that  is,  characteristic  boundaries). 

Problems  of  this  type  were  first  studied  for  linear  differential  equations 
in  two  dimensions  (that  is,  A(x,u)  = (a(x,y)  ,b(x,y))  and  h(x,u)  = 
h(x,y)  u + g (x,y)  in  (1.1))  by  O'Malley  in  a series  of  papers  [l6-l8]. 

More  recently,  Butuzov  [3>^]  has  considered  this  linear  equation  with  Dirichlet 
boundary  data  in  a rectangle  in  order  to  examine  the  interaction  of  the 
boundary  with  the  characteristic  curves  of  the  "semi-reduced"  equation 
y,(a(x,y)  ,b(x,y))*(ux  , u ) + h(x,y)  u + g(x,y)  = 0 . Finally  van  Harten  in 
his  dissertation  [20]  and  in  [21]  has  discussed  certain  linear  and  nonlinear 
equations  of  the  form  (1.1)  in  dimensions  greater  than  two  when  the  parameter 
(j,  is  a power  of  e . 
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The  article  of  O'Malley  [18]  the  monograph  of  Lions  [15]  and  the 
dissertation  of  van  Harten  [20]  also  contain  interesting  discussions  of 


related  singular  perturbation  problems  as  well  as  many  additional  references 
to  the  literature.  In  particular,  O'Malley  treats  many  results  on  two- 
parameter  problems  for  ordinary  differential  equations  and  these  can  serve 
as  motivation  for  seme  of  the  theorems  given  below.  CUr  approach  here  is 
fashioned  for  the  most  part  by  our  previous  discussions  of  singularly 
perturbed  partial  differential  equations  of  the  form  (1.1),  (1.2)  for  the 

case  m,  = 1 which  are  contained  in  [12,13].  These  results  and  those  of 
the  present  paper  are  proved  using  the  comparison  techniques  of  Amann  [1] 

which  compel  us  to  restrict  the  right  hand  side  of  equation  (1.2)  to  be  at 
most  a quadratic  function  of  the  first  derivatives  of  u . 

Part  I 

Dirichlet  Problems 

2.  The  Problem  (7^).  We  consider  first  the  Dirichlet  problem 

e Au  = u A (x,u)  • vu  + h(x,u)  , x in  n , 

u(x  ; e , |JU ) = cp(x)  , x on  T = 9 n , 

where  0 cEW  (N  > 2)  is  a bounded  region  which  is  described  by  a function 
F : ]RN  -»  K1  in  the  sense  that  n = {x  : F(x)  < 0]  . Consequently  the  boundary 
r of  n is  simply  F-1(0)  . The  parameters  e and  p,  axe  assumed  to  be 
small  and  positive,  and  to  tend  to  zero  in  an  interrelated  fashion. 

We  will  study  the  behavior  of  solutions  of  (7^)  with  the  help  of  certain 
solutions  of  the  reduced  equation  h(x,u)  = 0 obtained  by  formally  setting 
e and  ^ equal  to  zero.  In  order  to  simplify  the  presentation  we  will  assume 
(without  great  loss  of  generality;  cf.  Remark  2.1  below)  that  h(x,0)  =0  in  Cl  . 
Thus  our  study  concerns  those  solutions  of  (7?-^)  which  are  close  to  zero  in  Cl  > 


. 
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but  which  change  rapidly  near  T in  order  to  satisfy  the  prescribed  boundary 
conditions.  This  rapid  change  is  termed  boundary  layer  behavior,  and  it  is  a 
characteristic  occurrence  in  man^  singular  perturbation  problems  (cf.,  for 
example,  [18]).  Since  we  are  looking  specifically  for  such  solutions  of 
(7^)  we  define  now  the  closed  domain  in  Q x JR1  in  which  these  solutions 
lie,  namely  fl  = n X {u  : |u|  <C  d(x)|  , where  d > 0 is  a smooth  function 
defined  in  Q such  that  | cp(x)  | < d(x)  < | cp(x)  | + 8 for  x in 
and  d(x)  < 6 for  x in  Tg  • (Here  8 > 0 is  a small  constant  and  we 
introduce  the  sets  r = [x  in  Q : dist  (x,  T)  < p}  and  T'  = [x  in  fi  : 
dist  (x,  r)  > p}  for  any  positive  constant  p .) 

In  terms  of  the  domain  £ we  can  now  state  the  smoothness  assumptions 
which  we  impose  on  the  functions  occurring  in  the  statement  of  the  problem. 
First  of  all,  each  of  the  components  a.  of  the  vector  function  A(x,u)  = 
(a1(x,u) , . . . , a^(x,u)  ) as  well  as  the  function  h are  assumed  to  be  of 
Holder  class  < a < 1)  with  respect  to  x and  of  class  with 

respect  to  u in  £ . Similarly,  the  functions  F and  cp  are  assumed  to 
be  of  Holder  class  C^2,Q!\]RN)  and  c^2,0!^(rg)  , respectively.  Finally  we 
assume  that  the  trivial  solution  of  the  equation  h(x,u)  = 0 is  stable  in 

2 

the  sense  that  there  exists  a positive  constant  m such  that  hu(x,u)  > m > 0 
for  (x,u)  in  £ ; cf.  [12,13]. 

We  are  now  ready  to  discuss  the  existence  of  solutions  of  (7J  ) and  their 
behavior  for  small  values  of  e and  y,  . Suppose  first  that  y is  asymptoti- 
cally smaller  than  ^/e  (y  < < ) , that  is,  y./\/e  0+  as  e -»  0+  . 

Then  we  expect  that  the  influence  of  the  term  y A ■ v u on  the  behavior  of 
solutions  will  be  negligible  compared  to  that  of  h . Indeed,  there  is  the 
following  result. 


0 £ 
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Theorem  2.1.  Under  the  above  smoothness  and  stability  assumptions  there 


exists  a positive  constant  eQ  such  that  the  problem  (^)  has  a solution 
u = u(x  ; e , y,)  of  class  C^2’°^(Q)  whenever  0 < e < eQ  and  0 < p,  <<  J e . 
In  addition,  for  x _in  fi  we  have  that 


I u(x  ; e , p.)  | < Kexp  [m^xJ/Ve  ] , 

where  K = max  | ©(x)  | and  0 < mn  < m L_1  with  L = max  ||  v F(x)  ||  . (Here 

r n 

N 

and  below  |(  • H denotes  the  usual  Euclidean  norm  in  B .) 

Proof.  To  prove  all  of  the  theorems  involving  Dirichlet  problems  we  will  use 

the  following  well  known  result  (cf.  for  example  [1]).  Namely,  for 

N 2 

3(x  , u , vu  , p,)  = p,  A(x,u)*  vu  + h(x,u)  or  p,  E b (x,u) u + h(x,u)  satisfying 

j=l  J J 

the  given  smoothness  conditions  in  jS  the  boundary  value  problem 
e fiu  = 3(x  ,u,vu,p,),x  in  fi  , u = cp(x)  , x on  r = & fi  , has  a solution 

u = u(x  ; e , p.)  in  £)  of  class  C^2’a\  fi  ) for  each  e > 0 and  p,  > 0 for 

which  there  exist  functions  u>  , co  of  class  Cv  ’ '(n)  such  that  co  < co  , 
co  | < cp(x)  < co  | , and  e A<o  > 3(x  ,co  , v<n  , p.)  , eAco  < 3(x  , co  , vco  , p,)  for 

r r 

x in  Cl  . Moreover,  this  solution  satisfies  co(x  ; e , pk)  < u(x  ; e , p)  .< 

co  (x  ; e , p.)  for  x in  fT  and  for  all  such  e > 0 and  p,  > 0 . (The  function 

od(co)  is  called  a lower  (upper)  solution  of  the  boundary  value  problem.)  Such 
inequality  techniques  have  been  used  by  Dorr , Par ter  and  Shampine  [7]  and  the 

author  [ll]  for  related  one-dimensional  problems,  and  by  Eckhaus  and  de  Jager 

[9],  Eckhaus  [8]  and  the  author  [12,13]  for  higher  dimensional  ones. 

Thus  with  regard  to  Theorem  2.1  we  define  for  x in  [J  the  functions 

co, co  by  co(x,e ) = -<o(x,e)  = -K  exp  [m^x)/,/  e ] , where  K and  are 

given  above.  Clearly  for  x on  T co(x  , c)  < cp(x)  < co(x  , e)  , and  it  is 
just  as  easy  to  see  that  the  differential  inequalities  are  satisfied.  For 
example,  consider  co  . (The  verification  that  co  satisfies  the  proper 
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inequality  is  similar  and  is  omitted.)  Then  we  have  that  (for  3(x,u,7u,|j)  = 
pA(x,u)  • vu  + h(x,u) ) 

eAo)  - 3(x,a),  va),p,)  = {-n^AF^/e  - m1  ||vF|| 

+ (\k/jT  ) m]L  (A  • v F ) + h(x,0) 

+ hu(x,§)}  | CD  | 

> {-m-J  AF  | ^77  -m2  L2 
- (g./yT  ) m1  I A • V F I ^ + m2  ] | cd  | 

>0  by  our  choice  of  m^ 

for  e and  g,  sufficiently  small,  say  0 < e < , since  \i/J t -*  0+  as 

e -»  0 . Here  (x,£)  in  £ is  the  appropriate  intermediate  point.  The 

conclusion  of  the  theorem  now  follows  from  the  differential  inequality 
theorem  just  quoted. 

Suppose  next  that  the  parameters  g and  J e are  asymptotically  comparable 
(p,~v/e  ) , that  is,  there  exist  positive  constants  l and  such  that 

0 <ll<  \i/JT  < l as  e -»  0+  . Then  the  term  g,  A • v u contributes  more 
to  the  description  of  the  behavior  of  solutions  of  (71, ) than  it  did  in 
the  previous  case.  The  precise  result  is  the  next  theorem. 


Theorem  2.2.  Under  the  above  smoothness  and  stability  assumptions  there 
exists  a positive  constant  eQ  such  that  the  problem  (7^)  has  a solution 
u = u(x  ; e , g)  of  class  C^2,0t^(n)  whenever  0 < e < and  g,  ~ J e.  . 

In  addition,  for  x _in  Q we  have  that 

| u(x  ; e , u ) | < K exp  [nig  F(x)/y  e ] , 

I- 
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yftere  “2  > 0 --  s-Hc-h  that,  ^ L2  + ®2-t  L 2 < m2  . Here  K and  L are  as 

before  and  L=  max  | A(x,u)  • v F(x)  I . 

& 

Proof..  Define  the  same  functions  <n,u>  as  in  the  proof  of  Theorem  2.1 
(with  mx  replaced  by  m2  ),  and  note  that 

eAu>  -3(x  ,o)  ,vc£,  p)  > { -n^  | AF  | - m^L2 

2 

- lm2  L^+  m } | <n  | 

> 0 by  our  choice  of 

for  e sufficiently  small,  say  0 < e < e 

— o 

Suppose  now  that  p is  asymptotically  larger  than  (p,  >>  J~) 

that  is,  J e /p  -*  0 as  (4,  -*  0 . Since  p is  much  larger  than  J~ 

we  expect  that  the  term  p A • v u will  determine  the  behavior  of  solutions 
of  (flx)  at  least  near  the  boundary  of  n • Consequently  we  must  examine 
the  behavior  of  the  characteristic  curves  of  the  "semi- reduced"  equation 
(8)  pA(x,u)  • Vu  + h(x,u)  = 0 . Three  distinct  cases  are  present ; cf.  for 
example  also  the  related  discussion  in  [6].  First,  all  of  the  characteristic 
curves  of  (g)  may  leave  Q along  T either  tangentially  or  nontangent ially,  that  is , 

(*)  A(x,u)  • 7 F(x)  > 0 for  (x,u)  in  = ifl  ( r&  x E1)  . 

(If  T is  itself  a characteristic  curve  then  A(x,u)  • vF(x)  so  on  r .) 

We  have  the  following  result. 

Theorem  2_^.  Under  the  above  smoothness  and  stability  assumptions  there  exists 
JLJBQ.sitive  constant  pQ  such  that  the  problem  (7^)  has  a solution 
u = u(x  * e > U ) whenever  0 < p < pQ  , /T  <<p  and  the  inequality  (*) 
obtains . In  addition,  for  x in  Q we  have  that 


u(x  ; e , ji)  | < Kexp  [n^FCx)/  J e ] + p(e,g,)  , 


where  K and  are  as  before,  and  p(e,p,)  is  a positive  function  of 

order  (p,/v/e  ) exp  [ -8 /J  e ] <<  1 . 

Proof.  Define  for  x in  Q the  functions  by  cn(x  ; e , p)  = 

-o>(x  ; e , p.)  = -(Kexp  [n^Ffr)//®-  ] + p(e  , p,))  , where  p(e,u)  = 

_2 

m px(e  , p,)  for  P;L(e  , p,)  = 

{(g,/y7" ) m | A • VF  | +m  | A F | J~  + m2  L2  } exp  [ -m  S/VT"  ] • 

X oo  x oo  X X 

Then  proceeding  as  in  the  proof  of  Theorem  2.1  we  have  that 


2 2 

e Au>  - 3 (x  > “ > V<n  , p.)  > { | AF  1^  L 

+ (p ,/Jt  ) m^^  (A  • v F)  + m2}  | cn  + p | 

+ P1(e  , p,)  • 

Now  for  x in  Tg  A*  vF  > 0 and  so  we  have  the  desired  inequality  since 
m^  < m L 1 . But  for  x in  Tg  the  term  {•  } | <£  + p | is  such  that 
| {•}  | | + P | < p^(e  , p,)  and  so  the  desired  inequality  follows  here  also. 

Suppose  next  that  all  of  the  characteristic  curves  of  (g)  leave  n 
along  f nontangentially,  that  is,  there  exists  a positive  constant  k such (**) 


(**)  A(x,u)- vF(x)  > k||  V F(x)  ||  for  (x,u)  in  . 


Then  the  conclusion  of  Theorem  2.3  can  be  strengthened  as  follows. 

Theorem  2X.  Make  the  same  assumptions  as  in  Theorem  2.3  with  the  condition 
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Note  that  since  ,/e  <<  p , e/p,  <<  V/T  , and  so  the  boundary  layer  estimate 
of  this  theorem  is  sharper  than  that  of  Theorem  2.3- 

Proof.  Define  for  x in  n the  functions  co  , a)  by  co(x  ; e , p)  = 

—2 

-oj(x  ; e , p)  = - (K  exp  [kF(x)  p/e]  + p(e  , p))  , where  p(e  , p)  = m"  p-j/e  , p) 

for  P;L(e  , p)  = [p  k I AF  | ^ + (p2/e)(k2  ||  v F 11^+.  k | A * y F } Kexp  [ -k  5 p/e] 

It  follows  that 

eAo>-  ;?(x  ,cd  ,vo>,  p)  > (-k  | 

+ (p^/e  ) ( -k2  ||  v F ||2  + k A • V F) 

+ m } | u>  + p | + p1(e  ,p)  . 

2 

Now  for  x in  Tg  A " v F > k ||  v F ||  and  so  the  desired  inequality  is  valid 
for  p sufficiently  small,  say  0 < p < pQ  . And  for  x in  Tg  the  term 
{•}  |cu  + p | is  bounded  in  absolute  value  by  ; consequently,  we  have  the 
desired  inequality  for  all  x in  fi  . 

Suppose  finally  that  neither  of  the  two  previous  cases  obtains,  that  is, 
A(x,u)  • v F(x)  ^ 0 for  (x,u)  in  . Geometrically  this  means  that  some 
or  all  of  the  characteristic  curves  of  (§)  are  entering  Q along  portions 
of  r • Roughly  speaking,  in  this  case  we  have  no  control  over  the  boundary 
layer  behavior  of  solutions  along  portions  of  T where  A-vF  is  negative. 
This  is  reflected  in  the  relatively  crude  estimate  in  the  following  theorem. 


Theorem  2.5.  Under  the  above  smoothness  and  stability  assumptions  there 


exists  a positive  constant  pQ  such  that  the  problem  (7^)  has  a solution 
u = u(x  ; e , p)  of  class  C^2,Q^(n)  whenever  0 < p < pQ  and  ^/e  <<  p . 

In  addition,  if  A(x,u)  ’vF(x)  ^ 0 in  then  for  x _in  n we  have  the 

estimate 
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I u(x  ; e , p)  | < K exp  [m  F(x)  /p]  , 

2 

where  m > 0 is  such  that  m < m /i^  with  1^  = max  | A(x,u)  • 7 F(x)  | 
i * & 
and  K as  above. 

Proof.  Define  for  x in  Cl  the  functions  cn  , cu  by  cn(x,p)  = -<c(x,p)  = 

- K exp  [ m^  F(x)/p  ] . Then  we  have  that 

eA<u  - 3 (x  ,<u  , 7<a  , p)  > {-m  | AF  | e/p, 

— j <*> 

-m32  ||  7 F ||2  e/p2 
2 

-m^L^  + m } | (n  | 

> 0 by  our  choice  of  m^ 

if  p is  sufficiently  small,  say  0 < p < pQ  , since  e <<  <<  p . 

It  is  possible  to  improve  the  estimate  of  Theorem  2.5  by  considering 
more  closely  how  the  characteristics  of  (§)  interact  with  the  boundary  of 
Q • Before  doing  this  in  the  next  section  we  make  several  remarks. 

Remark  2.1.  If  u = uq(x)  / 0 is  a smooth  solution  of  h(x,u)  = 0 

2 

satisfying  hu(x,u)  > m >0  for  (x,u)  in  j0(uQ)  = n X [u  : | u - uq(x)  | < d^x)] 

(where  d^  > 0 is  a smooth  function  such  that  | cp(x)  - uq(x)  | < d^(x) 

< | cp(x)  -uo(x)  | + 6 for  x in  r&/2  6111(1  d]_(x)  — 8 for  x in  ) 

then  the  above  theory  applies  with  minor  modifications.  This  can  be  seen  by 

considering  the  problem  (77, ) for  the  function  z = u - u (x)  . 

1 o 

Similarly,  if  the  operator  eA  is  replaced  by  e L where  L is  the 
general  linear  second-order  uniformly  elliptic  operator  in  N variables  then 
our  results  for  this  more  general  problem  are  qualitatively  the  same  as  those 
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Remark  2.2.  The  stability  assumption  that  (f)  h (x,u)  > m >0  for 
(x,u)  in  & is  only  one  of  several  related  assumptions  which  often  have 
to  be  made  in  o^.  .er  to  study  the  problems  in  this  paper.  The  reader  can 
consult  [12]  for  precise  statements  of  these  assumptions.  With  regard  to  the 
condition  (f)  we  note  that  it  is  sometimes  necessary  in  order  to  prove 
results  like  Theorems  2.1- 2.3  and  2.5  to  make  the  weaker  assumption  that 

(i)  hu(x,0)  > m >0  in  Q ; and 

T| 

(ii)  cp(x)  J h(x,s)ds  > 0 for  x on  r and 

0 

0 < 7)  < cp(x)  (or  cp(x)  < T1  < 0 ) . 

Such  restrictions  allow  larger  "boundary  layer  jumps"  | cp(x)  | ; cf . , for 
example,  [10],  [ 12 ] or  Example  (E21)  below. 

Remark  2.3.  Some  of  the  results  of  this  section  and  several  extensions 
are  already  known.  The  case  of  a two-dimensional  linear  Dirichlet  problem 

(7^)  , that  is,  A(x,u)  = (a(x,y)  ,b(x,y))  and  h(x,u)  = h(x,y)  u + g(x,y)  , 

was  studied  in  [l6,l8],  [3, it-]  and  [20].  In  addition,  the  nonlinear  problem 

(7lj)  with  n the  N - dimensional  unit  ball  and  ^ = Je  was  considered 

recently  in  [21].  These  papers  treat  the  problem  of  Constructing  the  complete 
asymptotic  expansion  of  the  solution,  a problem  which  we  have  not  addressed  here. 

Remark  2 A.  If  A(x,u)  ’vF(x)  > 0 for  (x,u)  in  then  a compactness 

argument  shows  the  existence  of  a positive  constant  k such  that 

p 

(**)  A(x,u)  • vF(x)  > k||  vF(x)  ||  in  • Similarly,  if  there  is  a point 
xQ  on  T such  that  vF(xq)  = 0,  a(xq  , u)  f 0 and  A(x,u)  • vF(x)  > 0 
(for  all  u of  interest)  in  a punctured  neighborhood  G of  xo  then  the 
inequality  (**)  obtains  in  G U [xq]  . Thus,  with  regard  to  Theorem  2.h, 
the  boundary  of  n is  allowed  to  have  finitely  many  corners,  that  is,  points 


at  which  vF  vanishes. 
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Remark  2.5.  The  following  idea  may  have  occurred  to  the  reader  in  the  case 
that  y/  e <<  pi  • Namely  suppose  one  first  solves  the  "semi-reduced"  problem 
consisting  of  the  equation  (§)  together  with  a boundary  condition  along 
a portion,  of  T and  then  uses  this  solution  to  approximate  the  solution 
of  (7^)  . Specifically  the  semi-reduced  problem  is 


(S6>) 


p,  A(x  , u)  • vu  + h(x,u)  = 0 , x in  fi  , 
u(x>M,)  = cp(x)  , x on  r^r  , 


where  is  the  portion  of  T on  which  A(x,u)  *vF(x)  < 0 for  all  u 

of  interest.  A sufficient  condition  for  the  solvability  of  fS&)  is  that 
h(x,0)  2 0 and  hu(x,u)  >m  >0  for  (x,u)  in  ; cf.  [l4j  for  a proof 
in  the  case  N = 2 . Then  the  solution  u = uQ(x,p,)  of  (§#)  satisfies  in 

n 

uQ(x,n)  = 6(1^  exp  [mF1(x)/M,] ) , 

where  K.  = max  | cp(x)  j and  F = F | . Consequently,  along  the  portion 

T]_  ri 

of  T where  A(x,u)  *vF(x)  <0  an  easy  estimation  shows  that  the  solution 
u of  (7J^)  satisfies  | u-  uq  = o(l)  as  ^ 0+  , a result  which  agrees 

with  the  estimate  of  Theorem  2.5 . On  the  other  hand,  along  the  portion  of 
T where  A(x,u)  • vF(x)  > 0 the  same  kind  of  argument  used  in  the  proofs 
of  Theorems  2. 1-2. 3 shows  that  | u-uq|  = &(K  exp  [F(x)/s/~  ])  . But  for 
such  x , uq  = o(l)  as  p,  -»  0 and  so  this  estimate  agrees  with  our  previous 
one.  All  of  this  shows  that  we  can  obtain  the  results  of  this  section  in  the 
case  that  ,/ e <<  y,  by  using  instead  a solution  of  the  semi-reduced  problem 
(§&)  . We  find  it  more  convenient  to  use  the  solution  of  the  "totally"  re- 
duced problem  h(x,u)  = 0 . 
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3 • Some  Extensions . We  consider  in  this  section  some  sharper  results 

than  Theorem  2.5  in  the  case  that  J e <<  p,  and  A(x,u)  *vF(x)  ^ 0 

along  certain  portions  of  r • Specifically  we  have  in  mind  the  situation 

where  a characteristic  curve  of  the  semi-reduced  equation  (s)  is  tangent 

to  an  isolated  noncorner  point  of  T . 

Suppose  then  that  xq  on  T is  such  a point,  that  is,  yFfic^)^  0 , 

A(xq  , u)  f 0 but  A(xq  , u)  • v F(xq)  = 0 for  all  u of  interest.  Let 

y(x,u)  = A(x,u)  'vF(x)  in  a neighborhood  G of  xq  defined  by 

G = [x  : g(x)  < 0}  . In  terms  of  y and  G , we  assume  that  v(x0  j u)  = 0 

and  that  y changes  sign  in  passing  through  xq  , that  is, 

y(x,u)  < 0 in  G ~\{xo3  where  G"  = G fl  [x  : A(x,u)  • v g < 0 } and 
+ + 

y(x,u)  > 0 in  G \{x^}  where  G = G fl  {x  : A(x,u)  • V g > 0 ) . 

Proceeding  as  in  the  previous  section  we  now  examine  the  solutions  of 
(?}  ) in  the  subsets  G^  of  G . For  simplicity  we  assume  that  the  point 
of  tangency  is  xq  = 0 . 

Lemma  3.1.  Suppose  that  u = u(x  ; e , p,)  is  the  solution  of  the  problem 

■f* 

(7j^)  in  the  domain  £1  . Then  in  the  subset  G _of  G where  y > 0 we 
have  that 

u(x  ; e , g,)  = &(KX  exp  [ kF(x)/v] ) , 

where  K-.  = max  | cd(x)  | and  v = J e for  those  x such  that 
1 + 

G nr 

Y = While  v = ©(e/(y^.))  for  those  x such  that  y >>./T/  • 

Here  k > 0 is  a known  constant  depending  on  y . 

Note  that  for  y = 0 (y  - ©(l)  , y > 0)  this  estimate  agrees  with  the  one  of 
Theorem  2.3  (Theorem  2.4)  but  it  provides  much  more  information  than  these 


results. 
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I 

Proof.  To  prove  the  lemma  it  is  sufficient  to  note  that  for  x in 

-J- 

G (In  the  function  cd(x,v)  = ^ exp  [ kF(x)/v]  with  k and  v appro- 
priately chosen  is  a barrier  function  for  (7^)  . That  is,  o>  is  an  upper 
solution  and  -oo  is  a lower  solution  of  (7^)  ; cf.  [5]  or  [9].  This 
follows  in  the  case  of  o>  = -u>  , for  example,  because 

eAu>-  3(x,o>,Vcu  , y,)  > [ -k  |AF  | e/v 

-k2  ||  v F ||2  e/v2 

2 

+ky(x,a))  y,/v  + m }a> 

> 0 for  k appropriately  chosen 

provided  v = for  all  x such  that  y(x,cd)  = T /y,  ) and 

v = S(e/(v|i))  for  all  x such  that  ,/T/y  <<  y(x,co)  . A similar  argument 
shows  that  tu  is  an  upper  solution,  and  so  | u(  x ; e , y,)  | < <jd(x,v)  in 

G+nn  . 

Consider  finally  the  subset  G of  G where  y < 0 . 

Lemma  3-2.  Suppose  that  u = u(x  ; e , y,)  is  the  solution  of  (7^)  in  the 
domain  . Then  in  the  subset  G~  _of  G where  y < 0 we  have  that 

u(x  ; e , y,)  = exp  [ kF(x)/v] ) , 

where  and  k are  as  before  and  v = for  those  x such  that 

| Y | = ©(^/T/y,)  while  v = ©(  | y | y.)  for  those  x such  that 

| V | >>  ©(A/ M>)  • 
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Again  note  that  for  | y | sufficiently  small  ( | y | = 6(1)  , Y < 0 ) 
this  estimate  agrees  with  that  of  Theorem  2.1  (Theorem  2.5).  However,  it 
provides  more  precise  information  than  these  results. 

Proof.  The  proof  follows  that  of  Lemma  3.1  once  we  note  that  the  "y-  term" 
in  the  corresponding  inequalities  is  harmful  when  | y(x  , ±o>  )|  >>  6 (v^"/m>  ) 
and  so  we  must  choose  v = ©(  | Y I ^ ) for  such  x . 

Before  considering  another  class  of  Dirichlet  problems  in  the  next  section 
we  make  some  remarks  about  these  lemmas. 

Remark  3.1.  The  estimates  of  Lemmas  3.1  and  3.2  give  the  following  qualitative 

picture  of  the  behavior  of  the  solution  near  the  point  of  tangency  xq  = 0 . 

Namely,  for  y(x,±a>)~±  ||x||  , at  the  outer  edge  of  G " the  boundary 

layer  is  quite  thick  (of  order  p ) but  decreases  like  ©(  ||  x ||  p ) as 

||  x ||  decreases  until  ||x||  = ©(./7/p,  ) after  which  the  boundary  layer 

thickness  is  of  order  ,/7  . The  boundary  layer  thickness  remains  of  order 

y7  after  x passes  through  xq  = 0 into  G until  ||x||>>  ,/T / p 

after  which  the  thickness  decreases  to  ©(e  / 1|  x ||  p)  as  ||x||  increases. 

Finally  the  thickness  shrinks  to  order  e/p  as  y > 0 increases  to  order  one. 

The  linear  equation  e (u  + u ) = pu  + u , for  (x,y)  in  the  unit 

xx  yy  y 

disk  centered  at  (1,0)  , serves  as  a simple  illustration  of  this  remark. 

Here  y(x,y,u)  = y and  the  points  of  tangency  are  (0,0)  and  (2,0)  . 

Remark  3.2.  In  Lemmas  3.1  and  3.2  we  studied  the  solution  of  (7|  ) in 

the  domain  & . This  is  justified  by  ovn  stability  assumption  which  guarantees 

the  uniqueness  of  the  solution  of  (7J^)  in  the  domain  f)  * cf.  [19]  • 

Remark  3.3.  The  ideas  of  this  section  can  be  extended  to  cover  cases  in 
which  a portion  of  T actually  coincides  with  a characteristic  curve  of 
(S)  and/or  in  which  f has  corners  which  interact  with  the  characteristics. 


Both  situations  would  arise  if,  for  example,  n were  a rectangle  in  K 


and  the  characteristics  were  horizontal  or  vertical  lines.  We  discuss  such 
a problem  in  Section  7« 

Remark  3.4.  We  note  finally  that  if  y < 0 (in  G ) and  p,  = 1 then 
Lemma  3.2  shows  that  the  only  way  to  have  a close  approximation  to  the 
solution  of  (7^)  in  G ~ is  to  make  zero.  That  is,  we  must  choose 

a solution  u = uq(x)  of  the  reduced  equation 

(c^)  A(x,u)  • v u + h(x,u)  = 0 

which  satisfies  the  prescribed  boundary  data  along  that  portion  of  T . It 

is  this  solution  that  we  use  to  approximate  the  solution  of  (7^)  and  to 

which  we  would  add  a boundary  layer  corrector  term  along  the  part  of  T 

where  u does  not  necessarily  satisfy  the  boundary  data  ; cf.  for  example 
o 

[9],  [l8j  Sec.  7]  and  [15].  This  remark  also  makes  sense  geometrically 
because  along  that  part  of  T where  y < 0 the  characteristic  curves  of 
(S^)  are  incoming  and  so  it  is  possible,  in  general,  to  find  a solution 
of  (g1)  assuming  the  prescribed  boundary  values  there. 

k.  The  Problem  (7^)  . We  consider  next  the  following  problem  whose  right 
hand  side  contains  a quadratic  nonlinearity,  namely 
N 2 

eAu  = n £ b . (x,u)  u + h(x,u)  , x in  0 , 

3=1  J x3 

u(x  ; e , y.)  = <p(x)  , x on  T = d 0 . 


The  functions 
functions  a^ 
are  as  before. 


b . are  assumed  to  have  the  same  smoothness  properties  as  the 
considered  in  Section  2,  and  the  remaining  functions  and  sets 
We  will  assume  again  that  h(x,0)  = 0 in  n and  that  the 
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trivial  solution  is  stable  in  the  sense  that  there  exists  a positive 

2 

constant  m such  that  h (x,u)  > m >0  in  the  same  domain  £ . 

u — 

Our  discussion  of  the  asymptotic  behavior  of  solutions  of  ( 7[ 
parallels  that  of  Section  2 for  (7^)  ; however,  the  quadratic  nonlinearity 
introduces  some  distinctive  changes,  as  will  became  apparent  shortly. 

Suppose  first  that  y,  is  asymptotically  smaller  than  «(ii<<  e)  , 
that  is,  y,/e -»  0+  as  e -»  0+  . The  quadratic  dependence  of  the  right  hand 
side  on  the  first  derivatives  of  u prompts  us  to  compare  y,  and  e instead 
of  y,  and  ye  ; cf . the  following  proof  of  Theorem  4.1.  Then  the  analog  of 
Theorem  2.1  is  the  following  result. 

Theorem  4.1.  Under  the  above  smoothness  and  stability  assumptions  there 

exists  a positive  constant  eQ  such  that  the  problem  (TJg)  has  a solution 

(2  q)  

u = u(x;e  , y,)  of  class  C'  ’ ' (o)  whenever  0 < e < eQ  and  0 < y < < e . 

' In  addition,  for  x in  n we  have  that 

| u(x  ; e , y)  ( < K exp  [m^xJ/yT  ] , 

where  K = max  | cp(x)  | and  0 < m < m L 1 with  L = max  ||  v F(x)  ||  . 

r 1 n 

Proof.  Simply  proceed  as  in  the  proof  of  Theorem  2.1  with  3(x,u, v u,y)= 

N 2 - 

y £ b.(x,u)u  + h(x,u)  and  o>  , cjo  as  defined  there.  (Note  that 

d-i  0 xj 

No  + 

p.  I b (x,oj)oj“  = &((yi/e  ) exp  [2m  F(x)/y7  ] ) = ©(p/e)  -»  0 as  e -*  0 

d-i  0 J 

1 

for  cu  = 03  , tn  . ) 

Next,  if  y,  is  asymptotically  equal  to  e (pi  ~ e)  , that  is,  if  there 
exist  positive  constants  JL  and  such  that  0 < < p/e  < t , then 

we  expect  that  a result  like  Theorem  4.1  will  obtain  with  m^  replaced  by 
another  suitable  constant.  The  precise  result  is  the  next  theorem  whose 
proof  mimics  that  of  Theorem  2.2  and  is  therefore  omitted. 
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Theorem  4.2.  Under  the  above  smoothness  and  stability  assumptions  there 
exists  a positive  constant  eQ  such  that  the  problem  (7^)  has  a solution 
u = (x;e,p)  of  class  C^2,0L\ci)  whenever  0 < e < cq  and  0 < < p/e  < l 

In  addition,  for  x in  Q we  have  that 

| u(x  ; e , p)|  < K exp  [n^  F(x)/,/7  ] , 

where  0 < < m (L2  + l K 1^)  with  K and  L as  above.,_and 

N 2 , 

L = max  | £ b (x,u) F (x)  | . 

1 & J-l  J J 

Suppose  finally  that  p,  is  asymptotically  larger  than  e(p>>  e),  that 

is  - -*  0+  as  p -♦  0+  . Then  the  effect  of  the  quadratic  terms  becomes 
* U 

really  noticeable  for  the  first  time,  and  we  must  exercise  some  care.  As 
with  the  analogous  situation  in  our  discussion  of  (7?^)  three  distinct 
cases  present  themselves.  Suppose  first  that 
N 2 

(*)  <p(x)  £ b (x,u)  F (x)  > 0 for  (x,u)  in  fl&  . 

j-1  J J 

Geometrically  this  means  that  the  solution  (7?g)  must  be  convex  (cp  > 0) 

N 2 

or  concave  (cp  < 0)  near  T , since  for  such  x Au~(p/e)  £ b (x,u)Fx(x) 

“ j=l 0 0 

has  the  sign  of  cp(x)  as  a result  of  (*)  • 

Theorem  4.3.  Under  the  above  smoothness  and  stability  assumptions  there 
exists  a positive  constant  such  that  the  problem  (7^)  has  a solution_ 

u = u(x  ; e , p)  of  class  C^2,a\n)  whenever  0 < p,  < , 0<e<<  p and 

the  inequality  (*)  obtains.  In  addition,  for  x in  n we  have  that 

- K exp  [m1  F(x>//7  ] < u(x  ; e , p ) < 0 (if  cp  < 0) 
or 

0<  u(x  ; e , p)  < K exp  [m^  F(x)//7  ] (if  cp  > 0)  , 
where  K and  m^  are  as  before. 
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Proof.  Suppose  for  example  that  cp  < 0 and  define  for  x in  fi  u>  5 0 

and  o>(x,e)  = - K exp  [ra1  F(x)//e"  ] . Clearly  co  is  an  upper  solution  of 

(7^)  and  it  follows  just  as  easily  that  co  is  a lower  solution.  Indeed 

N 2 

we  have  that  (for  3(x,u,  vu,y,)  = y,  Z b.(x,u)u  + h(x,u)  ) 

j=l  J Xj 

2 2 

e&O)  - 3 (x,u>,  >f-m1|AF|v/e  - m.  L 

+ m2  (y/e ) ( Z b . (x,u>)  F 2 ) o>  + m2  } | cn  | 

0=1  J j 

> 0 by  our  choice  of  m^ 

if  y,  is  sufficiently  small,  say  0 < y,  < y,Q  . This  follows  because  for 

N 2 , 
x in  rc  -Z  b.(x,o>)F  (x)  >0  by  assumption,  while  for  x in  F~  | cn  | 

6 j=l  J “ Xj  " & “ 

is  transcendentally  small,  that  is,  | u>  | = ©(e11)  for  all  n > 1 . 

Thus  the  conclusion  of  Theorem  4.3  in  the  case  that  cp  < 0 follows 

from  the  differential  inequality  theorem  quoted  in  the  proof  of  Theorem  2.1. 

The  case  cp  > 0 is  handled  similarly  with  o>  = 0 and 

o>(x,e)  = K exp  [ m^FjxJ/yT  ] . 

Suppose  now  that  the  inequality  (*)  holds  in  the  stronger  sense  that 
there  exists  a positive  constant  k such  that  for  (x,u)  in 

No  2 

Z b (x,u)  F (x)  < -k  ||  vF(x)  ||  (if  cp  < 0)  or 

d-i  3 A 

(**) 

N 0 0 

Z b.(x,u)  F_“  (x)  > k ||  vF(x)  ||  (if  <p  _>  0)  . 

j=l  J Xj 

We  expect  that  the  estimate  of  Theorem  4.3  can  be  sharpened  accordingly. 

This  is  the  content  of  the  next  theorem. 


I 
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Theorem  4.k.  Under  the  above  smoothness  and  stability  assumptions  there 
exists  a positive  constant  pQ  such  that  the  problem  (7^ ) has  a solution 
u = u(x  j e , p)  of  class  C^2,a\n)  whenever  0 < p < pQ  , 0<e  <<  p , 
the  inequality  (**)  obtains  and  ||vF(x)||>  0 for  x in  r&  . In  addition, 
for  x _in  n we  have  that 

-1  "klv"1 

v k K £n  [ -F  (x)  + e ] < u(x  ; e , p)  < 0 (if  cp  < 0) 


or 


-k,  v 


-1 


0 < u (x  ; e , p)  < - v k"1  K h\  [ -Fq(x)  + e 1 ] (if  cp  > 0)  , 

where  v = e/p  , F (x)  = F(x)/(max  | F(x)  | + 1)  and  0 < k < k K with 

n 

K as  above. 


Proof.  Suppose  again  that  cp  < 0 and  define  for  x in  n cn  h 0 and 


-k,  v 


-1 


-1  1 
cn(x,v)  = v k^  K In  G(x,v)  , where  G(x,v)  = -Fq(x)  + e . Clearly 

us  is  an  upper  solution  and  as  for  o>  , we  have  that 


eAcn  - 3 (x  ,cn  , yen  >p)>  -{  ! &F0  ve  G 


-1 


+ G"2  ||  vFo||2 
N 


+ k_1K  pv2(  E b (x,cn)F  2 ) G 

.i=l  J °,Xi 


-2 


+ m2  v -tnG]  k^  1 K 


N 


Now,  for  x in  Tg  -E  b^(x,oj)  Fq  x > k ||  v Fq  ||  by  assumption  and  so 

J ”1  J 

the  desired  inequality  follows  by  our  choice  of  k^  if  p is  sufficiently 

small,  say  0 < p _<  pQ  . Finally,  for  x in  this  inequality  also 

2 

obtains  because  the  positive  terra  -m  v in  G dominates  all  of  the  others. 
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Suppose  finally  that  the  functions  b.  , F and/or  cp  are  such  that 

J N 2 , 

the  inequality  (*)  does  not  obtain,  that  is,  cp(x)  £ b (x,u)  F (x)  0 

j=l  0 j 

in  £ • Then  as  in  the  analogous  situation  for  the  problem  (71  ) we  have 

6 x 

no  control  over  the  behavior  of  the  solution  of  (7 ^ ) near  T . The  most 
general  statement  we  can  make  about  this  case  is  contained  in  the  next 
theorem. 

Theorem  4.5.  Under  the  above  smoothness  and  stability  assumptions  there 
exists  a positive  constant  p^  such  that  the  problem  (7Jg)  has  a solution 

u = u(x  ; e , u)  of  class  C^2,0^(fi)  whenever  0 < p < pQ  and  0 < e <<  p • 

N ? - 

In  addition,  if  cp(x)  £ b.(x,u)F  (x)  £ 0 in  £R  , then  for  x in  Q we 

j=l  J 3 

have  that 

I u(x  ; e , p ) | < Kexp  In^Fix)/^  ] , 

where  0 < < m (KL^)"1^2  with  K and  as  above. 

Proof.  Define  for  x in  n the  functions  u>  , o>  by  cn(x,p)  = -cn(x,p)  = 

- K exp  [ m2F(x)/yp  ] . Then  in  the  case  of  o>  we  have  that 

eAu>  - 3 (x  ,<n  , 7 0) , p)  > ( -n^  | AF  1^  e/^p" 

-m2||  VF||2  e/p 

2 N 2 

+ nu  ( £ b (x,u>)  F )u> 

^ d-1  J J 

+ m2 } | u)  | 

>0  by  our  choice  of  im, 

if  p is  sufficiently  small,  say  0 < p < pQ  . Similarly  u>  satisfies  the 
opposite  inequality,  and  so  the  theorem  follows. 
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We  close  this  section  with  several  remarks. 


Remark  4.1.  It  is  possible  to  consider  a more  general  problem  than  (7^>)  , 
namely 

N 2 

eAu  = Ul  Z b (x  , u)  u + u A(x,u)  • v u + h(x,u)  , x in  fj, 
j=l  3 o 

u(x  ; e , ^ , u2)  = <p(x)  , x on  r , 

where  e , ^ and  U2  are  small  interrelated  parameters.  However  we  have 
studied  only  the  simpler  problem  (7^)  in  order  to  isolate  the  effects  of 
the  quadratic  nonlinearities  on  the  asymptotic  behavior  of  solutions.  Depend- 
ing on  the  relative  sizes  of  e , and  y,2  solutions  of  (7^)  behave 

roughly  like  the  solutions  of  (7^)  and  (7^)  studied  above. 

Remark  b.2.  In  order  to  prove  Theorem  4.H  we  had  to  make  the  assumption 
that  vF(x)  f 0 for  x in  rg  , that  is,  we  did  not  allow  the  boundary 
to  have  corners.  If  however  vF(xq)  = 0 at  an  isolated  point  xq  on  T 
then  the  conclusion  of  the  theorem  is  still  valid  provided  that  in  a neighbor- 
hood G of  x we  use  IC  exp  [ mF(x)//e  ] (for  K = max  1 cp(x)  | and 

° 1 Gnr 

m > 0 appropriately  chosen)  as  the  boundary  layer  term.  The  validity  of 
this  substitution  follows  essentially  from  the  proof  of  Theorem  U.3. 

Remark  U.3.  Often  the  estimate  of  Theorem  U.5  can  be  strengthened  in  a 
manner  analogous  to  the  strengthening  of  Theorem  2.5  in  Section  3.  Namely, 

in  a neighborhood  G of  a point  x^  on  T in  which  B(x,u)  = 

N 2 , . 

£ b.(x,u)F  (x)  changes  sign  for  all  u of  interest  we  can  give  sharp 

j=l  J XJ 

estimates  for  the  boundary  layer  term.  These  estimates  depend  of  course  on 
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the  sign  of  B relative  to  the  sign  of  the  boundary  data  cp  , and  they 
are  the  analogs  of  the  estimates  in  Lemmas  3-1  and  3 >2.  We  leave  their 
precise  formulation  to  the  reader. 

Part  II 

Robin  Problems 

5.  The  Problem  (7?^).  We  turn  now  t°  a discussion  of  results  for  the  Robin 


problem 


(fl3 ) 


e A u = p,  A(x,u)  • 7 u + h(x,u)  , x in  n - {x  : F(x)  < 0]  , 
P(x)u  + vF(x)  • 7 u = cp(x)  , x on  T = d 0 , 


which  are  similar  to  those  for  the  corresponding  Dirichlet  problem  (7^)  • 

The  function  p is  assumed  to  be  nonnegative  on  r and  of  class 

C(2,a)(r&)  while  the  functions  A , h , F and  cp  are  as  above.  The  only 

exception  is  that  we  assume  here  and  in  the  next  section  that  ||  7F(x)  ||  = 1 

for  x on  T = F-1(0)  , that  is,  in  particular,  T is  not  allowed  to  have  corners. 

Thus  the  exterior  unit  normal  at  a point  x on  Tj  n(x)  , is  equal  to 

7 F(x)/||  vF(x)  ||  = V F(x)  , and  so  7F(x)-7u=  n(x)-7u=|^  is  the  outward 

normal  derivative.  In  view  of  this  the  boundary  condition  in  ( 7[ can  be 

du 

written  in  the  more  familiar  form  pu  + - cp  . 

We  assume  as  usual  that  h(x,0)  =0  in  Cl  ; however,  the  non- Dirichlet 
boundary  condition  prompts  us  to  seek  solutions  of  (7^)  which  are  uniformly 
close  to  zero  in  fj  (cf.  for  example£l3J ) • Thus  we  will  look  for  solutions 
of  (7 ) in  the  domain  e = 0 X { u : | u | < 6 } in  n x B1  and  we  ask 
further  that  the  zero  solution  be  stable  in  the  sense  that  there  exist  a 


2 

positive  constant  m such  that  hu(x,u)  > m >0  for  (x,u)  in  g • Note 
that  since  5 > 0 is  arbitrarily  small  (but  independent  of  e and  u ) this 
stability  assumption  is  essentially  that  hu(x,0)  >0  in  0 ; this  condition 
is  much  less  restrictive  than  the  corresponding  condition  for  the  Dirichlet 
problems  (71  ) and  (7^)  . 

Our  results  for  the  problem  (7^)  are  so  similar  to  those  for  the 

problem  (7^)  that  we  will  only  prove  the  first  of  the  two  theorems  stated 

below.  Moreover,  for  ease  of  presentation  we  will  state  these  results  in 

abbreviated  form,  it  being  understood  that  the  domains  & and  are  replaced 

by  the  domains  g and  gg  respectively  and  that  the  definition  of  stability 

is  that  given  in  this  section.  (Here  ge  = gH(rc  xK1)  .) 

o o 

We  consider  first  the  situation  in  which  the  ^ - term  does  not  really 
affect  the  behavior  of  solutions. 

Theorem  5«1-  The  conclusions  of  Theorem  2.1,  (Theorem  2.2)  and  Theorem 

2.3  are  valid  for  the  problem  (7^)  with  K replaced  by  ^/T  Km^1(v/e  Kn^-1)  . 

Proof.  The  theorems  for  Robin  problems  will  be  proved  by  using  a variation 
of  the  differential  inequality  theorem  quoted  in  Section  2 and  due  to 

Amann  [l].  Namely,  for  jj(x  , u , v u , |j, ) = u A(x,u)  • V u + h(x,u)  or 
N 2 

E b.(x,u)u  + h(x,u)  satisfying  the  given  smoothness  conditions  in  g 

d-i  J i 

the  boundary  value  problem  eAu  =3(x  ,u,7u,|i)  , x in  0 , 

Pu  + |^  = cp(x)  , x on  r = an  , with  P > 0 , has  a solution  u = u(x  ; e , p,) 
in  £ of  class  c^,Q^(o)  £°r  each  e > 0 and  w > 0 for  which  there  exist 
boundary-related  functions  u>,a>  of  class  C^’°^(n)  such  that  cn  < o>  , 

Pu)  + a 0)/  an|^,_<cp<PcD  + a cn/an  , and  e A to  > 3(x,  cu  , 7«)  > (x)  5 e A ox  < 

3(x,  ox  , 705  , p.)  for  x in  n . Moreover,  this  solution  satisfies  cd(x  ; e , y,)  < 
u(x  ; e , n)  < co(x  ; e , |±)  for  x in  n . The  bounding  functions  co  , cd  are  said 


to  be  boundary-related,  if  there  exists  a function  u - u(x  ; e , p)  of 


class  C^2,°^(fi)  such  that  for  x on  T P(x)u 
x in  Ci  cd(x  j e , p ) < u(x  ;e,p)<cD(xje,p)  . 


x on  T P(x)u  + 1^-  = cp(x)  and  for 


Turning  now  to  the  proof  of  the  theorem  consider  just  the  case 
u«/7  , as  the  other  two  cases  follow  similarly.  Define  for  x in 
n the  functions  cd , cd  by  cd(x,s)  = -CD(x,e)  = Km^  exp  [m^  F(x)/v/T  ] . 

Clearly  cd  < cd  and  P cd  + V F • v cd  | r < cp<PaJ  + vF‘va)  | • (Recall 

that  p > 0 . ) In  addition,  we  have  that  (for  3(x  , u , v u , p)  = 
p A(  x,u)  • v u + h(x,u)) 

eAcB-  3 (x  ,cd  , V<d  , p)  > { - | AF^e  - n^/e  ||  V F || 

+ pA-vF  + m2m11  JT  ] ( m1  //e" ) | cd  | 

> 0 by  our  choice  of  m^ 

if  e is  sufficiently  small,  say  0 < e < eQ  • Similarly  we  can  show  that 
bAcd  < 3(x  ,<d  , vcd,  p ) in  n , and  so  the  conclusion  of  the  theorem  will 
follow  once  we  establish  that  our  functions  cd  , cd  are  also  boundary  - related. 
To  this  end,  let  $ be  a C^2 ’a^  - extension  of  cp  to  0 , that  is,  cp  |p  = cp  , 
and  define  the  function  u by  u(x,e)  = T (e  ) u (x)  ( exp  [ F(x)/t(g  ) ] - 1), 
where  T(e)  >0  is  a transcendentally  small  term.  Then  clearly  we  have 
that  cd<u<cd  in  Q and  Pu  + vF'Vu|r=p-0  + £|r=cp,andso 
indeed  cd,cd  are  boundary  - related.  (Recall  that  ||  V F ||  s 1 on  T .) 

If  now  the  p - term  is  sufficiently  large  and  either  the  inequality 
(**)  obtains  or  the  inequality  (*)  does  not  obtain  then  we  have  the 
following  result. 

Theorem  5.2.  The  conclusion  of  Theorem  2.k  (Theorem  2-5)  is  valid  for  the 
problem  (7L)  with  K replaced  by  (e/p)  Kk  1 (p  Km^  1 ) . 
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We  close  this  section  with  two  remarks. 


Remark  5-l«  The  estimate  in  the  conclusion  of  Theorem  5-2  for  the  case 
A <<  jj,  and  A(x,u)  • vF(x)  ^ 0 in  Gg  can  be  sharpened  in  the  same 
way  that  Lemmas  3*1  and  3.2  sharpen  the  estimate  in  Theorem  2.5-  Namely, 
it  is  not  difficult  to  see  that  the  conclusions  of  Lemmas  3*1  and  3-2  are 
valid  for  the  problem  (7^)  with  replaced  by  v k 1 . 

Remark  5-2.  The  Robin  problem  (7^)  can  be  studied  under  definitions  of 
stability  which  are  more  general  than  the  one  used  here.  The  interested 
reader  should  consult  [13]  where  a related  class  of  problems  is  treated. 


6.  The  Problem  (7J^)  . We  consider  finally  the  problem 


N 2 

eAu  = p,  2 b.(x,u)u  + h(x,u)  , x in  Q = [x  : F(x)  < 0}  , 


(\> 


j=l 


x . 
J 


P(x)u  + v F(x)  • v u = cp(x),  x on  r , 


where  all  of  the  functions  appearing  here  are  as  defined  in  Sections  4 and  5. 

The  results  for  (7J^ ) analogous  to  those  for  the  Dirichlet  problem  (7^) 

follow  directly  once  the  appropriate  changes  are  made  in  order  to  accommodate 

the  non- Dirichlet  boundary  conditions. 

We  assume  then  that  h(x,0)  s 0 in  Q and  that  there  exists  a positive 

2 

constant  m such  that  h (x,u)  > m >0  for  (x,u)  in  the  domain  £ defined 

u — 

in  Section  5«  We  will  only  prove  the  first  of  the  next  three  theorems  since 
the  arguments  for  the  other  two  are  quite  similar. 

The  first  theorem  deals  with  the  case  in  which  the  jj,  — term  does  not 
really  affect  the  behavior  of  solutions. 
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Theorem  6.1.  The  conclusions  of  Theorem  4.1,  ( Theorem  4.2)  and  Theorem 
4.3  are  valid  for  the  problem  (7^ ) with  K replaced  by 

J~e  Kn”1  (/7  Knig-1  ) . 

Note  that  domains  £ and  & of  Section  4 are  replaced  by  the  domains 
£ and  £s  respectively. 

Proof.  We  will  prove  this  result  only  in  the  case  that  p,  <<  since 

the  proofs  in  the  other  two  cases  are  very  similar.  Define  for  x in  Q 
the  functions  cn  , cn  by 

co(x,e  ) = - co(x,e  ) = - Km^  1 exp  [ m^F(x)/y/~e  J . 

Then  w < w , Pen  + v F * vco  lr  < cp  < pin  +vF  ' V co  | ^ , and  for 
N 2 

3(x,u,  vu  , p,)  = y,  E b.(x,u)u  + h(x,u)  we  have  that 
J=1  ° XJ 

e Acn  - ;j(x,  co  ,Vu>  , y,)  > £ - e | A F L2 

_ N P 

-m  (li/v/7)(  E b (x,co)F  ) | cn  | 

1 0=1  J 0 

+ y7  nigin^1  KnyV  e ) | m | 

>0  by  our  choice  of  m^ 

if  e is  sufficiently  small,  say  0 < e < eQ  . Similarly,  we  see  that 
eh(x>  -3(x,a),vcn,p,)<0  in  fi.  The  conclusion  of  the  theorem  now  follows 
from  Amann's  theorem  quoted  in  Section  5 once  we  observe  that  the  verification 
of  the  boundary  - relatedness  of  a)  , en  follows  exactly  as  in  the  proof  of 
Theorem  5.1. 

We  consider  next  the  case  in  which  <<  p,  and  the  strong  inequality 

(**)  of  Section  4 obtains  in  the  domain  . 
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Theorem  6.2.  Suppose  that  J e <<  g and  that  the  inequality  (**)  of 
Section  4 obtains  in  £„  . Then  the  conclusion  of  Theorem  4.4  is  valid  for 

— g — — - 

the  problem  ( [71 ^ ) with  the  In  - term  replaced  by 
v-tn  [l  ( K"1  - kF(x)  v'1  ) “1  } , 

where  l is  such  that  K > 1 and  v = e/g  • 

Suppose  finally  that  <<  g,  but  that  the  inequality  (*)  of 

Section  4 does  not  obtain  in  , then  there  is  the  following  result. 

N 2 

Theorem  6.3.  Suppose  that  Jz  < < g and  that  cp(x)  £ b . (x,u)  F (x)  ^ 0 

j=l  J *0  | 

for  (x,u)  _in  . Then  the  conclusion  of  Theorem  4. 3 is  valid  for  the 

problem  (7^)  with  the  exponential  term  replaced  by  g Km^  exp  [ m F(x)/g  ] . 

We  close  with  the  observation  that  it  is  possible  to  consider  Robin 

N 2 

problems  for  the  more  general  equation  eAu  = g £ b.(x,u)u  + 

0=1  J xj 

g^  A(x,u)  • vu  + h(x,u)  . Moreover,  the  estimate  of  Theorem  6.3  can  be 
improved  by  closer  analysis  of  the  type  used  in  Section  3 (cf.  also  Remark 
5.1). 

Part  III 
Examples 

7.  Some  Examples.  In  this  section  we  discuss  several  examples  which  illus- 
trate the  theory  developed  above. 

Example  7.1.  Consider  first  the  Dirichlet  problem 

e Au  = e A(x,u)  * v u + u , x in  fj={x:  F(x)  < 0}  , 

(El) 

u(x,6 ) = cp(x)  , x on  r = an  , 

| 

* i 

— - j 
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where  o is  a positive  constant  and.  the  functions  A , F and  cp  have  the 
usual  smoothness  properties.  This  problem  has  a solution  u = u(x,e)  which 
is  unique  by  the  maximum  principle  (cf.  [l^])  for  each  e > 0 sufficiently 
Rmal  1 since  the  functions  o>  , (jd  defined  by  co  = -co  = -K  for  K = max  | cp(x)| 

r 

are  bounding  functions.  Moreover,  for  x in  Q we  have  the  (rather  crude) 
estimate  | u(x,e ) | < K . In  order  to  obtain  sharper  estimates  we  apply  the 
theory  of  Section  2. 

N 

Suppose  first  that  the  region  0 is  the  unit  ball  in  B centered  at 
0 , that  is,  F(x)  = \ ( II  x ||2  - 1)  . If  a>  j then  the  specific  form  of 
A is  immaterial  and  Theorem  2.1  or  2.2  shows  that  the  solution  of  (El) 
satisfies  in  C3 

(7.1)  u(x,e ) = <&(  K exp  [ F(x)/y"e  ) ) . 

However,  if  a < ^ then  we  know  that  it  is  the  interaction  of  A with  T 

which  determines  the  width  of  the  boundary  layer  along  T . For  example,  if 
2 2 2 

A(x,u)  = u x then  A(x,u)  * vF(x)  = u |jx||  > 0 , and  so  Theorem  2.3  im- 

plies that  the  solution  of  (El)  also  satisfies  the  estimate  (7.1)  for  this 
range  of  a . Similarly  if  A(x,u)  = kx  for  a positive  constant  k then 

O 

A(x,u)  'vF(x)  = k ||  x ||  and  we  deduce  from  Theorem  2.k  that  the  solution  u 
satisfies  the  sharper  estimate  in  fi 

(7.2)  u(x,e)  = ©(  K exp  [ kF(x)/e  1-a  ] ) . 

Finally  if  A(x,u)  is  such  that  A(x,u)  * x ^ 0 in  fig  then  the  most  general 
estimate  we  can  give  for  u is  that  of  Theorem  2.5,  namely  for  x in  n 


u(x,e ) = ®(  Kexp  [ F(x)/e°  ] ) . 
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For  such  an  A it  may  be  possible  to  give  sharper  estimates  like  those  in 
Lemmas  3.1  and  3.2  along  the  portions  of  T where  A(x,u)  * x > 0 or 
A(x,u)  • x > 0 . 

N 

Suppose  finally  that  Q is  the  unit  cube  in  E centered  at  0 , 
Np 

that  is,  F(x)  = - n (l-x.)  for  | x.  | <1  . Then  T = dn  = F (0)  has 

N ^ ^ N 

2 ’ vertices  or  corners  {§r]  and  vF(§  ) = 0 for  1 < r < 2 . Never- 


theless if,  for  example,  A(x,u)  = u x then  A(x,u)  ’vF(x)  - 

*N 


2 N 2 N 2 N- 1 2 . 

u x • 2(x  n (i-x.)  , x n (i-x  ),...,  x n (i-x  ))  = 
X J ^ J=1  J 1N  -5=1  J 


j=2  o - j-x  o " j=l 

p p N p p N ^ ' 2 2 N-l  o 

= 2u  {x  n (i-x.)  + xf  n (i-x  ) + ...+  x n (i-x  )}>o 
1 j=2  . J 2 j=l  J j=i  J 

jfe 


and  so  for  all  a > 0 the  estimate  (7.1)  is  valid.  While  if  A(x,u)  - x 

1 12  — 
and  a < 2 then  A(x,u)  * vF(x)  > ||  VF||  in  n and  so  the  estimate 

(7.2)  is  valid  with  k = ^ . 


Example  7 .2.  Consider  now  the  Robin  problem 


e Au  = eCT  A(x,u)  'vu  + u^-u,  x in  n , 


(E2) 


U + In  = <p(x)  » x on  r » 


,N 


where  a is  again  a positive  constant  and  Q is  the  unit  ball  in  E 

— = x * y u is  the  outward  normal  derivative 


centered  at  0 ; consequently. 


dn 


of  u . Suppose  that  the  function  A satisfies  our  smoothness  assumptions. 

3 

We  note  first  that  the  reduced  equation  h(u)  = u - u = 0 has  the  three 
solutions  0,1  and  -1  ; of  these,  ±1  axe  stable  while  0 is  unstable 
in  the  sense  that  h' (0)  = - 1 < 0 . The  theory  of  Section  5 applies  to  the 
stable  zero  1 (-1)  if  we  make  the  change  of  variable  u = u - 1 (u  = u + l), 
and  so  we  will  simply  quote  a theorem  from  this  section  with  the  tacit  under- 
standing that  the  appropriate  change  of  variable  has  been  made. 


30 


jf  0 > -|  or  if  0<  \ and  A(x,u)  • V F(x)  - A(x,u)  • x > 0 in 
e8  = [u  : | u-  1 | 5 &)  then  '^aeorein  5-1  implies  that  the  problem  (E2) 

has  a solution  u = u+(x,e)  such  that  in  Q 

(7. 3+)  u+(x,e ) = 1 + &(,/  e/2  K exp  F(x)/ v/"e  ] ) 

for  K+  = max  | q>(x)  - 1 | . Similarly,  for  these  choices  of  o and  for 
A(x,u)-x>ro  in  fig  = r&x  [u  : | u + 1 | < 5}  if  o < £ this  theorem 
implies  the  existence  of  a second  solution  u = u_(*,«)  such  that  in  Q 


(7*3  ) u_(x,e)  = - 1 + &(y  3/2  K exp  [v/r2F(x)/v/e  ] ) 

for  K = max  | cp(x)  + 1 1 • 

Supposernext  that  o < \ . If  A is  such  that  there  exists  a positive 
constant  k for  which  A(x,u)  * x > k ||  x ||2  in  4 and/or  H then  Theorem 
5.2  implies  that  the  problem  (E2)  has  a solution  u+  and/or  a solution  u_ 

such  that  in  ft 

0A±)  U (*,«)  =±1  + © (h  _1  fi1_0  exP  & k P(x)/e^”°  ) ) • 

Finally,  If  A(x,u)  • x i o in  e+8  and/or  t~  then  there  are  still  two 

solutions  u+  and  n_  of  (E2)  s hovever,  the  best  general  estate  »e  cm 

give  is 

u (x,«)  = ±1+ © (e°K±exp  [ F(x)/ea  J ) • 

A closer  analysis  may  allow  us  to  give  sharper  estimates  of  the  form  (7-3  ) 
or  (7.4*)  for  this  case  if  A(x,u)  * x > 0 along  portions  of  T (cf.  Remark 

5.1). 

Before  passing  to  another  example  we  note  that  a similar  equation  can 
be  used  to  illustrate  Remark  2.2.  Namely,  consider  the  Dirichlet  problem 
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e&u  = e°  A(x,u)  'Vu  + u-u^,  X in  O , 

(E2#) 

u(x,e)  = cp(x)  , x on  T . 

3 

Then  the  reduced  equation  h(u)  = u-u  =0  has  0 as  its  only  stable 

p 

solution  and  h'(u)  = l-3u  >0  provided  | u | < l/,/3  . Thus  for  all 
boundary  values  cp  such  that  max  | cp(x)  | < 1/V  3 we  can  apply  the  theory 

r 

of  Section  2.  However  this  restriction  on  cp  can  be  weakened  if  we  use 
the  integral  condition  (ii)  in  Remark  2.2,  that  is, 

J1  o 

cp(x)  J (s-sJ)ds  > 0 for  all  T]  in  (0,cp(x)]  or  [cp(x),0)  . 

0 

71  3 2 2 

For  instance,  if  cp(x)  > 0 then  J (s  - s"5)  ds  = -j-  (2  - T]  ) > 0 provided 
I _ 0 

0 < < J 2 , with  a similar  result  for  cp(x)  < 0 . Consequently  for  those 

cp  such  that  maoc  | cp(x)  | < J 2 we  can  actually  apply  the  theory  of  Section  2. 

r 

Example  7.3.  In  order  to  illustrate  the  theory  for  the  quadratic  problems 
(7^)  and  (7^)  we  consider  an  equation  which  models  a first-order  gaseous 
reaction  that  occurs  inside  of  a porous  catalyst  pellet,  namely  P -»  n Q , 
where  n is  the  stoichiometric  factor.  If  u is  the  normalized  concentration 
of  P and  0(~  n - 1)  is  the  volume  change  modulus  then  a possible  conservation 
law  governing  this  reaction  is  (cf.  [2]). 

(7-5)  V,[(l  + 0u)~1vu]=i2u  , x in  n . 

2 

Here  5 is  the  Thiele  modulus  which  measures  the  relative  importance  of 

reaction  versus  diffusion  and  n c ]R  is  the  idealized  region  in  which  the 

reaction  takes  place.  Suppose  further  that  the  reaction  produces  an  increase 

in  volume  (that  is,  0 > 0 ) and  that  the  reaction  is  diffusion  - limited 
2 

(that  is,  $ >>  1)  . Then  equation  (7-5)  can  be  cast  into  the  more  familiar 


form 


eAu  = e0(  1 + 0u)  1||  vu||2  + u(l  + 0u)  , 
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where  e = $ < < 1 . 

If  we  assume  that  a unit  concentration  of  P is  prescribed  on  the 
surface  T of  the  region  fi  then  the  problem  to  be  solved  is  one  of  the 
form  (7^)  , namely 

eAu  - e0(l  + 0 u)  -1  ||vu||2+  u(l  + 0u)  , x in  n , 


u(x,e ) s 1 j x on  T • 

For  g,  = e0  , (i  = 6(e)  <<  Jl.  since  0 = &(l)  , and  h"(u)  = u(l  + 0u) 
satisfies  h' (u)  > 1 for  u > 0 . (Since  u measures  concentration  we 
are  only  interested  in  nonnegative  solutions  of  (E3)  •)  Therefore  Theorem  4.1 
implies  that  the  problem  (E3)  has  a solution  u = u(x,e)  for  each  6 > 0 


such  that  in  Q 


(7.6) 


0 < u(x,e ) < exp  [ F(x)/y  e ] , where  f = F 


F_1(0)  . 


The  assumption  that  u h 1 on  T is  tantamount  to  the  assumption  that 
there  is  no  resistance  to  the  transfer  of  reactant  from  the  bulk  flow  to  the 
surface  of  the  pellet.  A more  realistic  assumption  is  that  there  is  in  fact 
a small  amount  of  resistance,  and  so  we  could  consider  the  more  general  problem 


(E  3 ' ) 


eAu=e0(l  + 0u)1||vu||2  + u(l  + 0u),x  in  n. 


u + ( gh)  _1  |^  = 1 , x on  T 


Here  gh  is  the  resistance  coefficient  which  is  known  as  the  Sherwood  number. 
This  time  Theorem  6.1  implies  that  the  problem  (E3  ')  has  a solution  u = u(x,e) 
for  each  e > 0 such  that  in  n 


(7.7)  0 < u(x,e  ) < ./e  gh  exp  [ F(x)/,/e  ] 
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Note  that  if  Sh  - ©(l/y/ e )>>  1 then  the  resistance  to  transfer  is 
negligible  and  the  estimate  (7-7)  reduces  to  {7.6). 

The  estimates  {7.6)  and  (7-7)  show  that  the  assumption  of  diffusion  - 
limitation  ($  >>  1 or  e <<  1)  implies  that  most  of  the  reaction  is 
confined  to  the  surface  of  the  pellet.  This  phenomenon  has  been  observed 
experimentally  (cf.  [2]). 

Example  7-^.  As  our  final  example  we  consider  the  linear  problem 
eAu  = y.ux  + u , (x,y)  in  Cl  = ((x,y ) : F(x,y)  < 0)  , 

(E4) 

u(x>y»e>M,)  = cp(x,y)  , (x,y)  on  r , 

where  is  the  square  of  length  2 centered  at  (0,0)  , that  is, 

E(x,y)  = -(l-x  )(l-y  ) . This  problem  for  ^/e  < < p,  was  studied  by 
O'Malley  [18  ; Sec.  7]  who  gave  estimates  for  the  solution  u in  the 
horizontal  strip  [(x,y)  :-l<x<l,-l<-l  + 5<y<l-6<l}.  He 
was  unable  however  to  analyze  u near  the  corners  of  Cl  and  along  the 
sides  y = ± 1 since  these  straight  lines  are  the  characteristic  curves 
of  the  semi-reduced  equation  ^ u^  + u = 0 . Using  our  theory  we  can  dis- 
cuss the  behavior  of  solutions  of  (E4)  throughout  n . 

First  of  all,  the  problem  (E4)  has  a unique  solution  for  all  e > 0 
and  p,  > 0 since  the  functions  cu  , ou  defined  by  cn  = -o>  = - K for 
K = max  | cp(x,y ) | are  bounding  functions  and  since  the  equation  satisfies 

r _ 

a maximum  principle.  If  p,<<yc  or  p,  ~ J e then  Theorems  2.1  and  2.2 
imply  that  this  solution  satisfies  the  estimate  in  n 


^(x,y;e,pi)  = ®(K  exp  [ F(x,y)/y7  ])  . 
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Suppose  however  that  ^/e  <<  ^ ; then  we  must  examine  the  function 
v(x,y,u)  = (1,0)  'VF  = 2x(l  - y ) . Since  y^O  for  (x,y)  in 
[ -1,1]  x [-1,1  ] the  most  general  estimate  in  Q we  can  give  for  u is 
(cf.  Theorem  2.5) 

u(x,y  ; e , = &(K  exp  [ F(x,y )/ p,  ] ) . 

Now  let  us  apply  the  results  of  Section  3*  First,  we  look  in  a neighbor- 
hood of  the  right  hand  boundary  x = 1 . For  all  y in  [-1+6,  1-6] 

(for  6 > 0 small  and  independent  of  e and  p)  a calculation  shows  that 

2 

there  is  a positive  constant  k = k(6)  such  that  y>  k ||vF||  , and  so 

Lemma  3.1  implies  that  for  such  x and  y 

u(x,y;e,p)  = &(K  exp  [ kF(x,y)  p/e  ] ) . 

This  estimate  agrees  with  that  of  O'Malley.  For  those  values  of  y in 
[ -1  , -1+  6)  and  (1-6,1]  we  still  have  that  y > 0 and  so  this  lemma 
implies  that  the  solution  u satisfies  there  the  estimate 

(7.8)  u(x,y  ; e , p ) = G>(  Kexp  [ F(x,y)/v] ) , 

where  v = -/e  for  y such  that  1-y  = ©(v/e /p)  and  v = ®(e/((l-y  ) p)  ) 

for  y such  that  1-y  >>  J e / ^ . Secondly,  we  look  in  a neighborhood 

of  the  left  hand  boundary  x = -1  . Here  y <0  and  so  Lemma  3*2  implies 
that  the  solution  u satisfies  the  estimate 


(7.9)  u(x,y  ; e , p,)  = ©(  Kexp  [ F(x,y)/V] ) 

where  v = J t for  y such  that  1-y  = ®(./e/n,)  and  v = ®((l-y  )p,) 
for  y such  that  1-y  >>  J e / |i  . Note  that  for  y in  [-1+6,  1-6] 
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the  estimate  (7*8)  becomes 

u(x,y;e,m)  = &(  Kexp  [ F(x,y)/y,  ] ) ; 

this  result  agrees  with  that  of  O'Malley. 

It  remains  to  consider  tne  upper  and  lower  boundaries  of  n . We  can 
restrict  our  discussion  to  the  upper  boundary  y = 1 and  note  that  our 
results  apply  to  the  lower  boundary  by  reflection.  Let  us  look  then  in  a 
neighborhood  of  y = 1 . For  x > 0 the  function  y is  nonnegative 
and  so  the  estimate  (7.8)  applies  to  u for  such  x and  y ; similarly, 
for  x < 0 we  have  that  y < 0 and  so  the  estimate  (7*9)  is  valid  for 
these  x and  y . 

Finally  we  remark  that  in  a neighborhood  of  the  corner  points  ( -1  , ±1  ) 
(where  y < 0 ) the  estimates  (7-9)  hold,  while  in  a neighborhood  of  the 
corner  points  ( 1 , ±1  ) (where  y > 0 ) the  estimates  (7-8)  are  valid. 
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